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1 
MAGNETIC RESONANCE IMAGING OF 
VENTILATION AND PERFUSION IN THE 
LUNG 


CROSS-REFERENCE TO RELATED 
APPLICATION 


This patent document claims priority of U.S. Provisional 
Patent Application No. 61/420,554, filed Dec. 7, 2010, 
entitled “IMAGING VENTILATION AND PERFUSION 
IN THE LUNG USING MRI”. The entire content of the 
before-mentioned patent application is incorporated by ref- 
erence as part of the disclosure of this application. 


STATEMENT REGARDING FEDERALLY 
SPONSORED RESEARCH OR DEVELOPMENT 


This invention was made with government support under 
grants EB005970, HL080203, and HL081171 awarded by 
the National Institutes of Health (NIH), along with a 
National Aeronautics and Space Administration (NASA) 
grant NCC9-58 #PF02103, awarded through the National 
Space Biomedical Research Institute (NSBRI). The govern- 
ment has certain rights in the invention. 


BACKGROUND 


This patent document relates to magnetic resonance 
image (MRI) technologies. 

The lungs are primarily responsible for gas exchange, as 
well as a number of other functions. Ventilation of the lungs 
refers to the intake and removal of air in/out of the lungs, 
which involves the mixing of inspired gas (inhaled gas that 
is humidified and at body temperature) with alveolar gas for 
gas exchange of oxygen (O,) and carbon dioxide (CO,) in 
the alveoli of the lungs. Perfusion refers to the delivery of 
blood to capillary beds of a biological structure, e.g., alveoli 
in the lungs. 

MRI is a medical imaging technique that applies a mag- 
netic field and a pulse of radio frequency (RF) energy to 
produce an image used for imaging internal biological 
structures, e.g., the lungs. MRI is based on the property of 
nuclear magnetic resonance (NMR). NMR is a physical 
property in which the nuclei of atoms absorb and re-emit 
electromagnetic energy at a specific resonance frequency in 
the presence of a magnetic field. The absorption and reemis- 
sion of energy can be dependent on the strength of the 
magnetic field and the magnetic property of the atoms (e.g., 
atoms whose nuclei possesses magnetic spin). 


SUMMARY 


Techniques, systems, and devices are disclosed for imple- 
menting a quantitative, non-injectable contrast MRI tech- 
nique to measure spatial ventilation-perfusion matching. 

In one aspect of the disclosed technology, a method using 
magnetic resonance imaging to characterize ventilation and 
perfusion in a lung includes acquiring a magnetic resonance 
(MR) image of the lung that includes MR data in a voxel of 
the MR image, which includes obtaining a breathing fre- 
quency parameter, determining a water density value, a 
specific ventilation value, and a perfusion value in at least 
one voxel of the MR image based on the MR data, which the 
water density value is used to determine an air content value, 
and determining a ventilation-perfusion ratio value that is 
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the product of the specific ventilation value, the air content 
value, the inverse of the perfusion value, and the breathing 
frequency. 

Various implementations of the above aspect can include 
one or more of the following features. The method can 
include the MR image being a 2D image. The method can 
include the MR image being a 3D image. The method can 
include determining the ventilation-perfusion ratio value 
voxel-by-voxel. The method can include acquiring the MR 
image using oxygen as an inhaled contrast agent. The 
method can further include determining an alveolar venti- 
lation value that is the product of the specific ventilation 
value, the air content value, and the breathing frequency. 
The method can include producing a spatial distribution of 
ventilation, perfusion, and ventilation-perfusion using the 
ventilation-perfusion ratio value, the alveolar ventilation 
value, and the perfusion value. The method can include 
producing an image that shows the spatial distribution of 
ventilation, perfusion, and ventilation-perfusion for each 
voxel associated with the acquired MR image, and the 
spatial distribution of ventilation, perfusion, and ventilation- 
perfusion can be used as an indicator of a healthy or a 
diseased lung. The method can include determining venti- 
lation-perfusion matching in the lung based on the ventila- 
tion-perfusion ratio value, and ventilation-perfusion match- 
ing can be used as an indicator of a healthy or a diseased 
lung. 

In another aspect, an MRI system to characterize venti- 
lation and perfusion in a lung includes an MRI machine and 
a processing unit. The MRI machine acquires an MR image 
of the lung. The processing unit is configured to command 
the magnetic resonance imaging machine to acquire the MR 
image that includes MR data in each voxel of the MR image, 
process the MR data to determine a water density value, a 
specific ventilation value, and a perfusion value in at least 
one voxel of the MR image, in which the water density value 
is used to determine an air content value, and determine a 
ventilation-perfusion ratio value that is the product of the 
specific ventilation value, the air content value, the inverse 
of the perfusion value, and a breathing frequency. 

Various implementations of the above aspect can include 
one or more of the following features. The MRI system can 
use oxygen as an inhaled contrast agent to acquire the MR 
image of the lung. The processing unit can also be config- 
ured to determine an alveolar ventilation value that is the 
product of the specific ventilation value, the air content 
value, and the breathing frequency. The processing unit can 
also be configured to produce a spatial distribution of 
ventilation, perfusion, and ventilation-perfusion using the 
ventilation-perfusion ratio value, the alveolar ventilation 
value, and the perfusion value. The processing unit can also 
be configured to produce an image that shows the spatial 
distribution of ventilation, perfusion, and ventilation-perfu- 
sion for each voxel associated with the acquired MR image. 
The spatial distribution of ventilation, perfusion, and venti- 
lation-perfusion can be used as an indicator of a healthy or 
a diseased lung. The processing unit can also be configured 
to determine ventilation-perfusion matching in the lung 
based on the ventilation-perfusion ratio value, and ventila- 
tion-perfusion matching can be used as an indicator of a 
healthy or a diseased lung. 

In another aspect, a computer program product compris- 
ing a nonvolatile computer readable storage medium having 
instructions stored thereon includes code for acquiring an 
MR image of a lung that includes MR data in a voxel of the 
MR image, code for determining a water density value, a 
specific ventilation value, and a perfusion value in at least 
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one voxel of the MR image based on the MR data, in which 
the water density value is used to determine an air content 
value, and code for determining a ventilation-perfusion ratio 
value that is the product of the specific ventilation value, the 
air content value, the inverse of the perfusion value, and a 
breathing frequency value. 

Various implementations of the above aspect can include 
one or more of the following features. The code for acquir- 
ing the MRI image of the lung can include code using 
oxygen as an inhaled contrast agent. The nonvolatile com- 
puter readable storage medium can further include code for 
determining an alveolar ventilation value that is the product 
of the specific ventilation value, the air content value, and 
the breathing frequency. The code for determining the 
alveolar ventilation value can include code for producing a 
spatial distribution of ventilation, perfusion, and ventilation- 
perfusion using the ventilation-perfusion ratio value, the 
alveolar ventilation value, and the perfusion value. The code 
for determining the alveolar ventilation value can also 
include code for producing an image that shows the spatial 
distribution of ventilation, perfusion, and ventilation-perfu- 
sion for each voxel associated with the acquired MR image. 
The nonvolatile computer readable storage medium can 
further include code for using the spatial distribution of 
ventilation, perfusion, and ventilation-perfusion to indicate a 
healthy or a diseased lung. The nonvolatile computer read- 
able storage medium can further include code for determin- 
ing ventilation-perfusion matching in the lung based on the 
ventilation-perfusion ratio value. The code for determining 
ventilation-perfusion matching can further include code for 
using ventilation-perfusion matching to indicate a healthy or 
a diseased lung. 

In another aspect, a method implementing a proton mag- 
netic resonance imaging technique using oxygen as a con- 
trast agent to quantify regional specific ventilation in a lung 
with magnetic resonance imaging includes acquiring T1 
weighted images of a slice of a lung of a subject, in which 
the images are voluntarily respiratory gated by having the 
subject take a single breath following each image and then 
relax to functional residual capacity in time for a subsequent 
image, and detecting alterations in the images associated 
with a regional rise in oxygen concentration that reflect 
specific ventilation. 

Various implementations of the above aspect can include 
one or more of the following features. The method can 
include acquiring T1 weighted images without a hyperpo- 
larized gas. The T1 weighted images can be 2D images. The 
Tl weighted images can be 3D images. The acquiring T1 
weighted images can use oxygen as an inhaled contrast 
agent. The method can further include producing an image 
that shows specific ventilation for each voxel associated 
with the acquired MR image. The method can include the 
specific ventilation being used as an indicator of a pulmo- 
nary pathological condition. 

The disclosed technology can provide one or more of the 
following advantages. For example, the described methods, 
devices and systems can provide a fully quantitative non- 
injectable contrast agent proton MRI technique to measure 
spatial ventilation-perfusion (V_,/Q) matching and spatial 
distribution of ventilation and perfusion. The disclosed MRI 
techniques can be performed using an MRI system such as 
a standard 1.5 T clinical scanner, which can allow routine 
clinical imaging of patients. Implementation of the disclosed 
techniques can allow quantitative, noninvasive evaluation of 
important aspects of lung function. Implementation of the 
disclosed techniques can eliminate contrast injection and/or 
radiation hazard. This can be important for repeated evalu- 
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ations, such as for research applications, and for screening 
patients with genetic predispositions to lung disease and 
with diseases such as collagen vascular diseases that have 
pulmonary manifestations. 


BRIEF DESCRIPTION OF THE DRAWINGS 


FIG. 1 shows an exemplary MRI-based ventilation-per- 
fusion matching system. 

FIGS. 2A-2D show diagrams of an exemplary MRI-based 
ventilation-perfusion method to characterize spatial V_,/Q 
matching from MR images. 

FIG. 3A shows an exemplary quantitative image of 
regional specific ventilation. 

FIG. 3B shows an exemplary quantitative image of 
regional air content. 

FIG. 3C shows an exemplary quantitative map of regional 
alveolar ventilation. 

FIG. 4A shows an exemplary quantitative image of per- 
fusion. 

FIG. 4B shows an exemplary quantitative measure of 
regional ventilation-perfusion ratio. 

FIGS. 5A and 5B show exemplary plots of the spatial 
distribution of V,,, Q and V,,/Q. 

FIG. 6 shows an exemplary plot of specific ventilation of 
a lung unit. 

FIG. 7A shows an exemplary plot correlation delay time. 

FIG. 7B shows an exemplary plot of the relationship of 
correlation delay to specific ventilation. 

FIG. 8 shows an exemplary time series of signal intensity 
for a single voxel. 

FIG. 9A shows an exemplary correlation delay time map 
in a sagittal slice of the lung. 

FIG. 9B shows an exemplary corresponding specific 
ventilation map. 

FIG. 10 shows an exemplary graph of specific ventilation 
for each lung region. 

FIGS. 11A, 11B, and 11C show exemplary graphs of 
specific ventilation plotted versus the vertical distance from 
the most dependent portion of the lung. 

Like reference symbols and designations in the various 
drawings indicate like elements. 


DETAILED DESCRIPTION 


Techniques, devices and systems are described for imple- 
menting a quantitative, non-injectable contrast MRI tech- 
nique to measure spatial ventilation-perfusion matching. 

A primary function of the lung is gas exchange. The 
efficiency of gas exchange can be determined by matching 
of delivery of fresh gas by alveolar ventilation (V_,) to local 
perfusion (Q). V_.,/Q matching can be affected by passive 
mechanisms such as anatomical matching of the branching 
airway and blood vessel structures and by active mecha- 
nisms such as hypoxic pulmonary vasoconstriction. V_,/Q 
matching can be used to characterize processes in which 
regions of the lung that receive fresh gas also receive 
deoxygenated capillary blood. Therefore, precise V_,/Q 
matching within the millions of gas exchange units can be 
considered a fundamental mechanism to determine the effi- 
ciency of the gas exchange process, which can provide 
information on health and disease. 

The measurement of ventilation-perfusion (V_,/Q) match- 
ing can be used as an indicator to distinguish between 
healthy or diseased lungs. For example, in normal or healthy 
lungs, the matching of ventilation to perfusion is relatively 
uniform, e.g., the mean and variance of V,/Q ratios are 
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substantially identical for both ventilation and perfusion, 
and regions of high ventilation also receive high perfusion, 
with regions of low ventilation receiving low perfusion. For 
example, optimal V_,/Q matching in the lung can be asso- 
ciated with tight, narrow, clusters of matched V,, and Q or 
with broader, highly correlated distributions of V_, and Q. In 
some cases, a broadening of the distributions with a main- 
tained V_,/Q matching can be considered as ‘stressed’ V_,/Q 
matching—although matching is preserved, the lung may be 
at risk. For example, the lung can be considered overbuilt for 
many of the demands of daily life (except for the most 
strenuous exercise), and a region of loss of function may be 
compensated for by the remaining normal lung, which may 
delay diagnosis. Measurement of the separate spatial distri- 
bution of V,, and Q as well as overall V,/Q can allow 
detection of this compensation. Therefore, it can be impor- 
tant to measure the spatial distribution of V_, and Q as well 
as overall V_,/Q matching. 

V_,/Q mismatch can be associated with lung disease, such 
as airway predominant diseases including asthma and dis- 
eases affecting the pulmonary circulation including pulmo- 
nary arterial hypertension (PAH). Lung disease is a leading 
cause of morbidity and mortality in the world. Various types 
of lung disease can be characterized by disruption of V_,/Q 
matching, resulting in gas exchange inefficiency. For 
example, decreased V_,/Q matching can be a hallmark of 
chronic obstructive pulmonary disease (COPD), with areas 
of both high V,/Q ratio (deadspace) and low V,,/Q ratio 
described, and an overall decrease in V_,/Q matching. Simi- 
larly, pulmonary hypertension, asthma, pulmonary edema, 
pulmonary fibrosis and adult respiratory distress syndrome 
can also be characterized by V_,/Q mismatch. 

For example, in patients with advanced lung disease, a 
diagnosis can be clinically observable, which can be con- 
firmed with some diagnostic tools. Despite this, existing 
diagnostic tools such as spirometry can lack sensitivity and 
specificity, which may limit their utility for early diagnosis. 
For example, spirometry is not defined as abnormal until on 
average the patient has lost 20% of lung function. This can 
suggest that the normal compensatory mechanisms of the 
lung can prevent early detection of abnormalities, and can 
delay diagnosis and treatment. For example, a highly inva- 
sive and sparsely available diagnostic tool known as mul- 
tiple inert gas elimination technique (MIGET) has been 
shown to be capable of identifying markedly abnormal V_,/Q 
matching in patients with mild COPD (GOLD stage 1); 
while in the same patients, the changes in spirometry were 
shown to be mild and arterial oxygen and carbon dioxide 
levels were largely normal. This disruption of V_,/Q match- 
ing in these patients can suggest that a quantitative technique 
to measure V,/Q matching can offer the ability to detect 
pre-clinical disease in COPD. Similarly, patients with pul- 
monary hypertension often do not present for medical care 
until the disease is well advanced. Pulmonary hypertensive 
patients also may have marked disruption of V_,/Q matching 
at the time of diagnosis, suggesting that increased V_,/Q 
mismatch can be an early feature of the disease. The 
disclosed technology can include a quantitative measure- 
ment of spatial V_,/Q matching that can be used to enhance 
early disease identification and provide a noninvasive tool 
for monitoring disease progression and response to therapy, 
as well as quantitatively address research questions related 
to disease mechanisms. 

The disclosed technology described in this patent docu- 
ment can be used to directly measure whole lung and spatial 
V_,/Q matching in a manner that is fully quantitative, non- 
invasive, without radiation exposure, and clinically relevant 
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(e.g., applicable and accessible to human patients). Imple- 
mentation of the disclosed technology can provide the 
ability to measure spatial V_,/Q matching and spatial distri- 
bution of ventilation and perfusion, which can, for example, 
be performed on a standard 1.5 T clinical scanner at the level 
of the gas exchange unit. Measured V,/Q matching and 
spatial distribution of V, and Q can be used to evaluate 
mechanisms of lung disease, injury, response to treatment 
and repair. For example, the disclosed techniques may be 
used in screening populations of patients at risk of devel- 
oping lung disease and for monitoring response to treatment 
in those already afflicted. Routine clinical imaging of 
patients can be performed using the disclosed technology, 
e.g., allowing quantitative and noninvasive evaluations of 
the of lung function without concerns of radiation hazards 
and contrast injections. 

In one aspect of the disclosed technology, three quanti- 
tative non-injectable contrast proton MRI techniques, that 
when implemented together, can allow regional ventilation, 
perfusion and spatial V_,/Q matching to be directly mea- 
sured, e.g., in the entire lung in humans. The measured 
values can be based on regions of the lungs. For example, 
MR images of the lung of a subject can be partitioned into 
three gravitational regions, corresponding to thirds of the 
lung based on equal vertical extent, e.g., a dependent region 
(region which is lowest with respect to the direction of the 
gravity vector); an intermediate region (region between 
dependent and non-dependent regions); and the non-depen- 
dent region (region which is highest with respect to the 
direction of the gravity vector). One technique of the dis- 
closed technology can include a T2* multi-echo, fast gra- 
dient echo (mGRE) sequence for the rapid measurement of 
lung density (fractional water) in g (of water)/em>. Another 
technique of the disclosed technology can include a specific 
ventilation imaging (SVI) technique, which can use oxygen 
as a contrast agent for the measurement of regional specific 
ventilation, e.g., the dimensionless ratio of tidal volume 
(fresh gas) to functional residual capacity (residual gas at the 
end of a normal exhalation) of each lung element. Another 
technique of the disclosed technology can include an arterial 
spin labeling (ASL) technique, which can directly quantify 
regional pulmonary blood flow (perfusion) in units of mL (of 
blood)/min/em*. The disclosed MRI technology can apply 
appropriate scaling such that the measurements provided by 
the three exemplary techniques (e.g., MGRE, SVI, and ASL 
techniques) can be processed to provide quantitative spatial 
V_,/Q matching and spatial distribution of V_, and Q, e.g., on 
a voxel by voxel basis. For example, by combining the 
measured specific ventilation data with proton density mea- 
surements, alveolar ventilation (V_,) can be determined, such 
that when combined with perfusion measurements, can 
provide a fully quantitative measure of V_,/Q used to deter- 
mine V,/Q matching. For example, the disclosed MRI 
technology can measure T1 enhancement as a function of 
breath number to calculate an oxygen wash-in, and from 
that, be used to determine specific ventilation. Separate 
proton density measures can then be used to calculate V,, 
from specific ventilation, and combined with separate mea- 
sures of perfusion giving a fully quantitative measure of 
ventilation-perfusion matching. 

FIG. 1 shows an exemplary MRI-based ventilation-per- 
fusion matching system (100) for providing a quantitative, 
noninvasive measure of spatial V,/Q matching. FIG. 1 
shows one aspect of the exemplary system (100) that can 
include a magnetic resonance imaging (MRI) machine 
(110), which can be controlled by a MR image processing 
unit (120). 
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The exemplary MRI machine (110) can be used in the 
system (100) to implement a MRI-based ventilation-perfu- 
sion matching characterization process under the control of 
the exemplary MR image processing unit (120). MRI 
machine (110) can include various types of MRI systems, 
which can perform at least one of a multitude of MRI scans 
that can include, but are not limited to, T1-weighted MRI 
scans, T2-weighted MRI scans, T2*-weighted MRI scans, 
spin (proton (‘H)) density weighted MRI scans, diffusion 
tensor (DT) and diffusion weighted imaging (DWI) MRI 
scans, Tlp MRI scans, magnetization transfer (MT) MRI 
scans, real-time MRI, functional MRI (fMRI) and related 
techniques such as arterial spin labeling (ASL), among other 
MRI techniques. 

The exemplary MR image processing unit (120) can 
include a processor (121) that can be in communication with 
an input/output (1/O) unit (122), an output unit (123), and a 
memory unit (124). MR image processing unit (120) can be 
implemented as one of various data processing systems, 
such as a personal computer (PC), laptop, tablet, and mobile 
communication device. To support various functions of MR 
image processing unit (120), the exemplary processor (121) 
can be included to interface with and control operations of 
other components of MR image processing unit (120), such 
as the exemplary I/O unit (122), the exemplary output unit 
(123), and the exemplary memory unit (124). 

To support various functions of the MR image processing 
unit (120), memory unit (124) can store other information 
and data, such as instructions, software, values, images, and 
other data processed or referenced by processor (121). 
Various types of Random Access Memory (RAM) devices, 
Read Only Memory (ROM) devices, Flash Memory devices, 
and other suitable storage media can be used to implement 
storage functions of memory unit (124). The exemplary 
memory unit (124) can store MRI data and information, 
which can include subject MRI image data including spatial 
and spectral data. MRI machine system parameters, data 
processing parameters, and processed parameters and data 
that can be used in the implementation of a V,/Q matching 
characterization. Memory unit (124) can store data and 
information that can be used to implement a MRI-based 
V_,/Q matching characterization process and that can be 
generated from a MRI-based V_,/Q matching characteriza- 
tion algorithm and model. 

To support various functions of the MR image processing 
unit (120), the exemplary I/O unit (122) can be connected to 
an external interface, source of data storage, or display 
device. Various types of wired or wireless interfaces com- 
patible with typical data communication standards, such as 
Universal Serial Bus (USB), IEEE 1394 (FireWire), Blu- 
etooth, IEEE 802.111, Wireless Local Area Network 
(WLAN), Wireless Personal Area Network (WPAN), Wire- 
less Wide Area Network (WWAN), WiMAX, IEEE 802.16 
(Worldwide Interoperability for Microwave Access (Wi- 
MAX)), and parallel interfaces, can be used to implement 
T/O unit (122). I/O unit (122) can interface with an external 
interface, source of data storage, or display device to retrieve 
and transfer data and information that can be processed by 
the processor (121), stored in the memory unit (124), or 
exhibited on the output unit (123). 

To support various functions of the MR image processing 
unit (120), the output unit (123) can be used to exhibit data 
implemented by the exemplary MR image processing unit 
(120). The output unit (123) can include various types of 
display, speaker, or printing interfaces to implement the 
exemplary output unit (123). For example, the output unit 
(123) can include cathode ray tube (CRT), light emitting 
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diode (LED), or liquid crystal display (LCD) monitor or 
screen as a visual display to implement the output unit (123). 
In other examples, the output unit (123) can include toner, 
liquid inkjet, solid ink, dye sublimation, inkless (such as 
thermal or UV) printing apparatuses to implement the output 
unit (123); the output unit (123) can include various types of 
audio signal transducer apparatuses to implement the output 
unit (123). The output unit (123) can exhibit data and 
information, such as patient diagnostic data, MRI machine 
system information, partially processed MRI-based V,/Q 
matching characterization processing information, and com- 
pletely processed MRI-based V_,/Q matching characteriza- 
tion processing information. The output unit (123) can store 
data and information used to implement a MRI-based V_,/Q 
matching characterization process and from an implemented 
MRI-based V_,/Q matching characterization process. 

FIG. 2A shows a process flow diagram that describes 
implementation of an exemplary MRI-based ventilation- 
perfusion method to measure spatial V_,/Q matching from 
MR images, e.g., from the lung that can be used to assess 
lung function. Exemplary methods to characterize lung 
functionality using MRI techniques, such as water density 
measurement, specific ventilation measurement, alveolar 
ventilation measurement, perfusion measurement, ventila- 
tion-perfusion ratio measurement, and ventilation-perfusion 
matching, can be performed on systems and devices, such as 
the exemplary system (100) illustrated in FIG. 1. For 
example, FIG. 2A shows an exemplary method (200) to 
characterize spatial V_,/Q matching from MR image(s) that 
includes a process (210) to determine lung density (water 
fraction) values from an MR image of a lung in a desired 
target, a process (220) to determine regional specific venti- 
lation values from an MR image of the lung in the desired 
target, a process (230) to determine regional pulmonary 
perfusion values in an MR image of the lung in the desired 
target, and a process (240) to determine spatial ventilation- 
perfusion matching in the lung of the desired target, e.g., by 
combining and processing the resultant values process 
(210), process (220), and process (230). Process (240) can 
also include characterizing the spatial distribution of venti- 
lation and perfusion. Exemplary method (200) can imple- 
ment exemplary processes (210), (220), and (230) in the 
order shown in FIG. 2A, or in any other order not shown in 
the exemplary figure. 

Process (210) can include performing MR imaging on an 
MRI machine, e.g., MRI machine (110), using 2D or 3D 
multiple gradient-recalled echo (mGRE) sequence and/or 
other types of sequence. Process (210) can characterize 
regional lung water density (water fraction) by using an 
exemplary mGRE sequence technique to provide rapid 
acquisition of lung proton density data on a voxel by voxel 
basis. For example, lung proton density data can be acquired 
at rates ~8 seconds/slice. The described mGRE sequence 
technique can provide a measurement of regional lung water 
compensated for potential differences in T2*. 

The exemplary mGRE sequence technique can include a 
pulse sequence of a small flip angle slice selective excitation 
followed immediately by the collection of a single full line 
of Fourier data. The sequence can acquire data at 2 echo 
times, and the proton density of the lung can be determined 
by back extrapolating the signal to an echo time (TE), e.g., 
TE=0, by fitting data points to a single exponential. 

FIG. 2B shows a diagram of an exemplary 2D acquisition 
sequence to determine lung density (water fraction) values 
in process (210). In this example, process (210) can include 
process (211) that includes collecting 2D MRI images in 
sequential slices across the lung field, e.g., at 2 echo times. 
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Implementation of the exemplary 2D acquisition sequence 
on a slice with a thickness of substantially 15 mm can be 
covered over the entire lung in approximately 9 minutes 
(e.g., ~18 slices*10 sec per scan, 20 sec recovery). The 
exemplary 2D acquisition sequence to acquire MRI images 
in process (211) can include a single RF pulse followed by 
a single full line of Fourier data. It is noted that acquisition 
time is limited by breath-hold time (e.g., less than 12 
seconds in some patient populations. Process (212) includes 
phase encoding the acquired image(s) at a fixed echo time 
(TE) and repetition time (TR) until k-space is filled for a 
single image. The sequence can then repeat by filling 
k-space acquired at a second echo time while keeping TR 
constant. Process (213) includes processing the phase 
encoded data, e.g., by fitting the data points to a single 
exponential decay function, to determine lung water density, 
as well as T2*. The exemplary mGRE sequence technique 
exemplified in FIG. 2B can be adapted from a 2D acquisition 
sequence to a 3D acquisition sequence. For example, an 
exemplary 3D acquisition sequence can be implemented in 
which one 3D image is acquired at each echo time (e.g., two 
3D volumes), shortening acquisition time. For example, a 
single lung volume may be acquired in 6 seconds with an 
isotropic resolution of 8 mm* with the following sequence 
parameters: echo time=0.7 ms and 1.2 ms, field of view=30 
cm (A/P)x48 cm (R/L)x48 cm (S/D, repetition time=2.5 ms, 
flip angle=5°, matrix acquisition=40 (A/P)x64(R/L)x64 
(S/T), and can increase in signal-to-noise ratio (SNR). 

Process (220) can include performing proton MR imaging 
on an MRI machine, e.g., MRI machine (110). Process (220) 
can characterize specific ventilation (SV) and regional ven- 
tilation by using an exemplary specific ventilation imaging 
(SVI) technique. Specific ventilation (SV) is a measure of 
regional ventilation that can be determined by the ratio of 
fresh gas entering a lung region divided by its end-expira- 
tory volume. Specific ventilation (SV) can be defined as the 
ratio of the volume of fresh gas (AV) moving into a region 
of the lung to the end-expiratory volume (V,) of that region, 
SV=AV/V,. 

The exemplary SVI technique can be implemented in a 
single or multiple lung slice enabling two dimensional (2D) 
and three dimensional (3D) SV characterizations. The exem- 
plary SVI technique can use 100% oxygen as a contrast 
agent. The SVI technique can include oxygen delivery to 
and dissolution in lung tissue, which can shorten the longi- 
tudinal relaxation time (T,), increasing the local MRI signal. 
For example, T,-weighted images can be acquired while a 
subject breathes air and then 100% O, in an alternating 
sequence. Since the amount of regional dissolved oxygen 
can be associated with the amount of oxygen locally avail- 
able (e.g., determined by ventilation), the oxygen wash-in 
time for individual voxels can be determined. Voxels that 
have a high specific ventilation can reach the new equilib- 
rium faster than units that are less ventilated. A quantitative 
measure of the specific ventilation (e.g., local tidal volume/ 
functional residual capacity (FRC)) of that portion of the 
lung can be determined, for example, considering the time 
delay between the onset of the stimulus (breathing 100% 
oxygen) and the response for a particular voxel. 

FIG. 2C shows a diagram of an exemplary MRI imaging 
technique to determine to regional specific ventilation values 
in process (220). In this example, process (220) can include 
process (221) that includes acquiring 2D T, weighted MRI 
images, e.g., using an inversion recovery single shot fast 
spin echo sequence with images being acquired with a half 
Fourier acquisition. The exemplary process (221) can 
include using an inversion time of 1000 msec. Process (222) 
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includes producing a signal intensity time series of the 
acquired MRI images and determining the correlation delay 
time using a shifted cross correlation between the inspired 
fractional oxygen concentration and the acquired signal 
intensity time series. Process (223) includes processing the 
correlation delay time to determine specific ventilation val- 
ues, e.g., regional specific ventilation values, as described in 
further detail later in this patent document. Process (223) can 
also include producing a plot or map of the specific venti- 
lation values. The exemplary SVI technique is not limited to 
2D MRI acquisitions. The exemplary technique exemplified 
in FIG. 2C can be adapted from a 2D MRI technique to a 3D 
MRI technique. 

Exemplary techniques to implement 3D ventilation imag- 
ing can be accomplished in a reasonable acquisition time, 
e.g., by co-registering with the proton density images, 
specific ventilation can be converted to V,, as previously 
described. The exemplary SVI technique can map the 
change in signal intensity due to the local change in oxygen 
concentration effect on T,. This contrast change can be 
enhanced by adding a global inversion recovery pulse prior 
to acquisition. As a result, 3D sampling techniques can be 
implemented to acquire the entire lung field. For example, 
one approach includes improving the acquisition speed 
(short T, in the lung) by a single-shot GRASE acquisition 
which can allow multiple 2D images to be acquired follow- 
ing an inversion pulse. Another approach includes imple- 
menting highly constrained back-projection for time-re- 
solved MRI (HYPR) or vastly under-sampled isotropic 
projection reconstruction (VIPR) to significantly reduce the 
acquisition time and the echo time, thereby improving the 
overall SNR while allowing full 3D acquisition of the lung. 

Process (230) can include performing MR imaging on an 
MRI machine, e.g., MRI machine (110), using 2D or 3D 
calibrated arterial spin labeling (ASL) technique and/or 
other technique. Process (230) can characterize the distri- 
bution of pulmonary perfusion by using the exemplary ASL 
technique to provide analysis of how the normal distribution 
of perfusion can be disrupted by a variety of insults, such as 
hypoxia, exercise and interstitial pulmonary edema. The 
described ASL technique can be respiratory and cardiac 
gated and fully quantitative. 

The exemplary ASL technique can be used to quantify 
pulmonary perfusion using a 2D arterial spin labeling— 
flow-sensitive alternating inversion recovery with an extra 
radio frequency pulse (ASL-FAIRER) sequence with a 
half-Fourier acquisition single-shot turbo spin-echo 
(HASTE) imaging scheme. For example, two ECG-gated 
images (to collect data from a complete systolic ejection 
period) of each of the selected slices can be taken (e.g., 4 to 
6 seconds apart, depending on heart rate) during two suc- 
cessive apneas with the glottis open at the desired lung 
volume. The images can be preceded by tagging pulses that 
invert the magnetization of arterial blood outside the imag- 
ing plane in one case, and leave it undisturbed in the other. 
Pairs of images so obtained can be subtracted, leaving an 
image of pulmonary blood flow delivered during one cardiac 
cycle. The orientation of the imaged slices of lung can be 
varied as desired. One example includes the sagittal plane, 
and voxel sizes of approximately 1.6x1.6x15 mm (~36 mm? 
or 0.04 cm?) can be visualized. 

FIG. 2D shows a diagram of an exemplary MRI imaging 
technique to determine to regional perfusion values in 
process (230). In this example, process (230) can include 
process (231) that includes acquiring 2D MRI images, e.g., 
an ASL-FAIRER sequence with a HASTE imaging scheme 
in sequential slices. In FAIRER, an image is acquired by 
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applying a spatially selective inversion pulse, followed by a 
saturation pulse, followed by dephasing gradients. An 
adequate delay can allow relaxed spins to enter from adja- 
cent tissue and displace saturated spins. Imaging can begin 
after this delay. The imaging slice and inversion slice can be 
centered on the same plane to overlap. Complete inversion 
of the spins in the image plane can include making the width 
of the selective inversion pulse substantially wider than the 
image slice thickness. This can be repeated with the spatially 
selective inversion pulse replaced by a nonspatially-selec- 
tive inversion pulse. This can, for example, invert all of the 
spins in the volume covered by the RF body coil. The same 
delay can allow inverted spins to enter the imaging plane. 
Process (232) can include producing a perfusion-weighted 
image map, e.g., by subtraction of the spatially-selective 
tagged image from the non-spatially-selective tagged image. 
It is noted that the addition of the extra RF pulse and 
dephasing gradients, for example, can prevent the zeroing of 
signal intensity in certain voxels upon subtraction. Process 
(233) can determine regional pulmonary perfusion values by 
extracting perfusion data from the produced perfusion- 
weighted image map. 

An exemplary 3D ASL acquisition technique to acquire 
perfusion data in the lung can include collecting each 
tag/control pair in sequential slices such that the entire 
dataset encompasses the entire lung. Data for each sagittal 
slice can be collected during sequential breath holds at 
functional residual capacity moving from the medial chest 
adjacent to the mediastinum to the lateral chest wall. A 
complete 3D lung dataset may only take ~3 minutes to 
acquire with respiratory gating between repetition times, 
e.g., ~5 sec. This whole sequence can be repeated for signal 
averaging or to ensure image quality. A more complex, but 
more efficient, approach for 3D ASL data acquisition can 
include interleaving the ASL with the SVI acquisition, 
shortening the combined acquisition time for the two 
images. 

FIGS. 3A, 3B and 3C show exemplary quantitative 
images of regional specific ventilation (FIG. 3A) and 
regional air content (FIG. 3B) to produce a map of regional 
alveolar ventilation (FIG. 3C) in a single sagittal slice from 
the right lung of an exemplary human subject in the supine 
position. FIG. 3A shows quantitative image (301) that 
exemplifies regional specific ventilation in the exemplary 
slice by implementing process (220). FIG. 3B shows quan- 
titative image (302) that exemplifies regional air content in 
the exemplary slice by implementing process (210). FIG. 3C 
shows quantitative image (303) that exemplifies regional 
alveolar ventilation (V_,) determined on a voxel by voxel 
basis in the exemplary slice by implementing process (240). 
The exemplary color scale represents regional alveolar ven- 
tilation in L-min™'. The exemplary arrow (310) in FIGS. 3A, 
3B, and 3C indicates the gravitational vector. As seen in 
FIGS. 3A and 3C, the exemplary quantitative image values 
of both specific ventilation image (301) and regional venti- 
lation image (303) are greater in the gravitationally depen- 
dent region of the lung. This result may indicate, for 
example, that the nondependent region of the lung is already 
stretched and relatively poorly ventilated. 

FIGS. 4A and 4B show exemplary quantitative images of 
ASL measure of perfusion (FIG. 4A), which when combined 
with the measure of alveolar ventilation (303) can provide 
regional V_,/Q ratio (FIG. 4B). The exemplary data shown in 
FIGS. 4A and 4B were analyzed from the exemplary sagittal 
slice of the subject as shown in FIGS. 3A-3C. FIG. 4A 
shows quantitative image (401) that exemplifies ASL mea- 
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sured perfusion data in the exemplary slice by implementing 
process (230). The exemplary arrow (410) indicates the 
gravitational vector. 

FIG. 4B shows quantitative image (402) that exemplifies 
the regional ventilation-perfusion ratio in the exemplary 
slice by implementing process (240). In this exemplary 
implementation, the sagittal plane was chosen to demon- 
strate the well-known gravitational gradients in ventilation, 
perfusion and V_,/Q ratio. In the example shown in FIG. 4B, 
the regions of higher V_,/Q ratio are seen in the nondepen- 
dent region of the lung. In the example, regions consistent 
with large vessels can be mapped as regions of shunt V_,/Q 
(e.g. V_/Q ratio=0), as shown by arrow (420). The exem- 
plary arrow (410) indicates the gravitational vector. 

Process (240) can be used to determine spatial V_,/Q 
matching and spatial distribution of V_, and Q by implemen- 
tation of an exemplary MRI-based V_,/Q matching technique 
described herein. Specific ventilation can be measured in 
process (220), e.g., having units of mL (fresh gas)/mL 
(resident gas). Fractional lung density for any voxel can be 
measured in process (210), e.g., having units of mL (water)/ 
cm’. In this context, both blood and other intra/extra-cellular 
water can be considered as water. For example, a voxel in 
the lung can include a mixture of water (fractional density 1) 
and air (fractional density 0). Thus, the value of [1-density] 
can be a direct measure of the air content of a voxel, which 
can be measured in units of mL (resident gas)/cm*, e.g., 
within a small error resulting from any tissue not producing 
an MRI signal. 

The exemplary MRI-based V_,/Q matching technique of 
process (240) uses the combination of specific ventilation 
(e.g., measured using the SVI technique in process (220)) 
and density (e.g., measured at FRC using the mGRE 
sequence technique in process (210)) on a voxel by voxel 
basis to provide a measure of fresh gas input to that voxel, 
in units of [mL (fresh gas)/mL (resident gas)]:[mL (resident 
gas)/cm*], e.g., mL (fresh gas)/cm*. For example, quantita- 
tive data from an acquired SV image (e.g., image (301)) can 
be multiplied by the breathing frequency (f,) in units of 
min“', which can be imposed by the repetition time used in 
the SVI acquisition. The resultant value can be multiplied by 
the value of [1-density](e.g., quantitative data of air content) 
to provide a fully quantitative measure of alveolar ventila- 
tion (V,,), in units of mL (fresh gas)/cm*/min. V,, can be 
determined for each voxel, which is represented in Eq. (1): 


V4=SVIx(1-density)xfp (1) 


where SVI represents quantitative data from an acquired SV 
image; (1-density) represents quantitative data of air con- 
tent; and f, represents the breathing frequency. 

For each voxel, a V_,/Q ratio can be calculated by dividing 
the resultant V, image data (e.g., image (303)) by the 
quantitative value of perfusion in a voxel from an image 
(e.g., image (401)), e.g., from quantitative perfusion data 
measured using the ASL technique in process (230). V_,/Q 
can be determined for each voxel, which is represented in 


Eq. (2): 


V_,/Q=[SVIx(1-density)/ASL]xfp (2) 


where V_,/Q represents a quantitative ratio of alveolar ven- 
tilation to perfusion and ASL represents quantitative perfu- 
sion data. From Eq. (2), V_,/Q is represented in units of [mL 
(fresh gas)|/[mL (blood)]. 

Process (240) can further include a process to evaluate 
spatial distribution of V,, Q and V,,/Q ratios, the relation- 
ship of o(V_,/Q), the standard deviation (SD) of V_,/Q ratios 
in logarithmic space, and o(V,,) and o(Q), the SD of the 
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individual ventilation and perfusion distributions (e.g., in 
2D) logarithmic space. For example, after calculating the 
voxel by voxel V_,/Q ratio, some 5000-7000 points may be 
left, which can then be probed in a variety of ways. For 
example, plots of regional V_, versus Q can be generated and 
the histograms of the individual distributions (V_,, Q, and 
V_,/Q) can be quantified. The correlation (r) between 
regional V_,/Q and V,, and Q can be calculated. The variance 
in the V_,/Q heterogeneity can be related to the variance in 
V_, and Q as described by Eq. (3) and Eq. (4) in logarithmic 
form: 


0/0) =0(V 4P+0(Q)P-2 7-0 V.)(Q) (3) 
0° log(V4/0)=07 log V4+07 log 0-2-0 log Vo 


log Q (4) 


where ris the correlation coefficient between V_, and Q. Eqs. 
(3) and (4) show that different sets of o(V_,), o(Q), and r can 
produce substantially identical o(V_,/Q). 

FIGS. 5A and 5B show two exemplary plots of the spatial 
distribution of V,, Q and V,/Q for two healthy normal 
subjects. Plot (501) shown in FIG. 5A features spatial 
distribution of V,, Q and V,/Q of one normal healthy 
subject. Plot (502) shown in FIG. 5B features spatial distri- 
bution of V,, Q and V,/Q of the other normal healthy 
subject. The diagonal red line (510) represents a V_,/Q ratio 
of 1. Voxels above line (510) have a V_,/Q ratio above 1 (e.g., 
‘higher’ V_,/Q ratio). Voxel points below line (510) have a 
‘lower’ V_,/Q ratio. Plot (501) and plot (502) exhibit the 
overall V_,/Q distribution slightly above 1, in which the 
majority of points range between a V_,/Q ratio of 0.5 and 3. 
The range of V,/Q ratios can be considered important 
because it can be shown that this is the range of V_,/Q ratios 
consistent with normal blood oxygenation and thus repre- 
senting optimal V_,/Q matching. 

The exemplary method (200) can be implemented using 
system (100) or other MRI systems. In one example imple- 
menting the disclosed technology, MRI data can be collected 
from the exemplary subjects using a 1.5 Tesla Signa HDx 
TwinSpeed MRI system (e.g., General Electric Medical 
Systems, Milwaukee, Wis., USA). MR images can be taken 
from, for example, a single sagittal slice of the lung. A 
subject can be placed supine in the MRI scanner wearing a 
face mask (Hans Rudolph, Kans., USA), equipped with a 
non-rebreathing T-valve. The inspiratory port of the T-valve 
can be connected to a remote controlled valve allowing the 
operator to change from room air to gas contained in a large 
gas-tight bag filled from a gas cylinder in the control room. 
The inspiratory paths can be matched in airflow resistance to 
avoid any alteration in FRC when switching between them. 
The outlet of the T-valve can be connected to a long (e.g., 6 
m), low resistance expiratory line, e.g., leading out of the 
scanner room, where expired tidal volume can be measured 
(e.g., using a ParvoMedics Metabolic Measurement System, 
(Sandy, Utah, USA)). Reference phantoms can be placed on 
the anterior chest wall so as to be within the field of view of 
the scans, permitting absolute quantification of blood flow 
and lung density. The subject can lie on the posterior element 
of an 8-channel cardiac coil, and the anterior elements are 
placed directly on the chest wall. 

The exemplary MRI V,,/Q matching technique can be 
implemented in 2D and 3D. 3D implementations of the 
disclosed technology may include additional considerations, 
e.g., accurate image registration at a sub-voxel resolution. 
Other considerations that can exist in 3D implementations 
can include effects on image registration from motion and 
tissue deformation as lung volume at FRC has small varia- 
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tions in subjects. Additionally, large blood vessels may 
represent shunt V_,/Q ratios (e.g., an absence of ventilation 
that maps to a V,,/Q ratio of zero) if a voxel represents an 
imaged space entirely filled with blood. However, since 
these large vessels are conduits, this does not represent a true 
shunt, which is an important contribution to pulmonary gas 
exchange abnormalities. The disclosed MRI V_,/Q matching 
technique can include mechanisms to eliminate the contri- 
bution of large vessels from the image to accurately quantify 
true shunt. For large airways, these can be mapped as 
deadspace (regions that are ventilated but not perfused). 
Also, 3D implementations can involve longer data acquisi- 
tion times than 2D, for example. 

An important advantage of the disclosed techniques is that 
there is no “dose” associated with performing measure- 
ments, as there is no radiation or contrast injection to a target 
subject. Thus, the disclosed technology can allow repeated 
measurements to be performed on a subject under multiple 
conditions. For example, this can further benefit biomedical 
and physiological research design because each subject can 
act as their own control. 

In another aspect of the disclosed technology, specific 
ventilation can be measured using proton MRI by using O, 
as a contrast agent. Exemplary implementations of disclosed 
techniques of this aspect can demonstrate a mapped distri- 
bution of specific ventilation in the lung, e.g., including the 
lung’s gravitationally dependent gradient. Disclosed are 
exemplary devices, systems, and methods to quantify spe- 
cific ventilation (SV). 

Specific ventilation is a dimensionless quantity that pro- 
vides a measure of how efficiently a given lung region is 
ventilated. For example, regional specific ventilation can be 
an important metric from a physiological standpoint, as well 
as also for clinical diagnosis and assessment of different 
pulmonary pathological conditions, e.g., asthma, emphy- 
sema, cystic fibrosis and lung cancer. 

Oxygen is weakly paramagnetic, and when in solution in 
lung tissues can produce a measurable decrease in the 
longitudinal relaxation time T1 of tissues, increasing the 
signal intensity of an appropriately timed inversion recovery 
proton MR image. An exemplary technique can include 
using proton MRI technique in which inhaled oxygen (O,) 
can be used as a contrast agent to verify the presence or 
absence of ventilation. Additionally, the disclosed technol- 
ogy can obtain an O,-enhanced MRI signal that contains 
more information than just the presence or absence of 
ventilation. For example, the change in O, concentration in 
lung tissues (determining the local change in T1) depends on 
the rate of change of the O, alveolar concentration, which is 
a function of the regional specific ventilation. Thus, mea- 
surement of the regional rate of change of the MRI signal 
can allow the quantification of specific ventilation: Follow- 
ing a change in inspired fractional oxygen content (F,O,), 
units that have higher specific ventilation reach the new 
equilibrium faster than units that have a lower specific 
ventilation. Thus, the rapidity of the change in the MRI 
signal in a particular voxel or region following a change in 
inspired F,O, is a quantitative measure of the specific 
ventilation of that portion of the lung. 

In one example, following a sudden change in inspired 
fraction of O, (F,O,), the rate of change of the alveolar O, 
concentration is a function of the local specific ventilation, 
e.g., lung units with higher specific ventilation reach a new 
equilibrium faster than units with lower specific ventilation. 
The time it takes to reach a new equilibrium is a measure of 
the local specific ventilation. This time can be reflected in 
the MR images. For a series of inversion recovery images 
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acquired with appropriate scanning parameters, the signal 
intensity change observed following a change in inspired gas 
can be determined by the change in relaxation time T1, e.g., 
a change driven by the local amount of oxygen in solution, 
which is a function of the local O, partial pressure. The time 
course of the regional MR signal intensity resulting from a 
change in inspired O,-concentration can therefore reflect 
local specific ventilation. For example, the rapidity of the 
increase in the MRI signal intensity following the onset of 
O, inhalation is directly related to specific ventilation (as 
shown in FIG. 6, discussed later). 

Exemplary Lung Unit Model and Simulated Data 

An exemplary model of a lung unit is disclosed that can 
be used to translate the MRI signal rise time into a quanti- 
fication of specific ventilation. The exemplary model of a 
lung unit corresponds to the voxel size measured with the 
MRI. The variable of interest in the exemplary model is the 
rate of equilibration (e.g., the rise time). For example, dead 
space, water vapor and end tidal partial pressure of CO, 
(PCO,), factors that change the steady state equilibrium 
partial pressure of O, (P,O,), but remain largely constant 
during, can be ignored in many instances, as described 
herein. The exemplary lung unit model describes the tem- 
poral behavior of the O, concentration from the air breathing 
equilibrium value to the new 100% O, steady state. 

Different ventilation-perfusion (V_,/Q) ratios can result in 
different steady state oxygen concentrations for both 
inspired air (e.g., ~range 50-100 mmHg, for V_,/Q ratios 
between 0.001 and 10) and 100% inspired O, (e.g., ~550 to 
600 mmHg, for the same range of V,/Q ratios). These 
exemplary steady state equilibrium conditions can provide 
adequate (and essentially similar) contrast for all reasonable 
V_/Q ratios. Further, the exemplary lung unit model is 
independent of the equilibrium or steady state conditions, 
e.g., it can be based on the rate of equilibration. For example, 
two lung units presenting identical V_,/Q ratios and different 
specific ventilations are still distinguishable from each other 
by the different rate of equilibration, with the one presenting 
a higher SV having a faster turn-over, and thus equilibrating 
faster. 

For the exemplary lung unit model, assumptions can be 
applied. For example, it can be assumed that a voxel 
dimension (1.6x1.6 mm in the imaging plane, for a 15 mm 
thick slice, volume ~40 mm*) is small enough such that 
concentrations inside the unit at end expiration can be 
considered uniform, and therefore each voxel can be treated 
as a single ventilatory unit. Also, it can be assumed that the 
time interval between two consecutive acquired images (>5 
sec) is long enough so that equilibrium between oxygen in 
the gas-phase and in solution in blood and tissues is attained 
by end expiration. Using either continuous (trumpet model) 
or discrete (asymmetrical branch point) models of the aci- 
nus, these exemplary approaches have shown that at the 
scale of existing resolution, oxygen concentration can be 
constant within a single respiratory pathway. Exemplary 
simulations in asymmetrical branch point models suggest 
differences in O, concentration can persist among parallel 
respiratory units, yet within each unit the O, concentration 
gradient can be essentially abolished, consistent with the 
exemplary assumption. 

In the exemplary simulation, subjects can self gate their 
breathing to the 5 s interval between consecutive image 
acquisitions. Oxygen dissolves in tissues very rapidly com- 
pared to this 5 s interval between images. For example, for 
a normal, room air inspired breath, it takes ~0.25 s for 
oxygen to diffuse through a 0.5 um thick capillary wall and 
reach its equilibrium with hemoglobin in the pulmonary 
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capillary, a process that includes dissolving in tissue as well 
as other processes. The process of simple dissolution in 
tissue thus can be assumed to be complete within this same 
time frame, as in the exemplary assumption. 

FIG. 6 shows illustrations and plots of specific ventilation 
of an exemplary lung unit model. Illustration (610) and (611) 
illustrate specific ventilation in the lung unit, in which V, 
represents the end-expiratory volume of the unit and AV 
represents the volume increase that occurs during inspiration 
in the lung unit. Each voxel (of an MR image) can be 
simulated as a single ventilatory unit, with end expiratory 
volume V,, to which inspiration transiently adds AV during 
each breath (as illustrated in FIG. 6 by illustration (610) and 
illustration (611)). The initial concentration in the unit at end 
expiration is denoted C,. C,, (n=1, 2, 3, .. . ) denotes the 
concentration in the unit at end expiration after breath (n 
denotes the breath number). 

For example, let C,,,,,," denote the concentration of oxy- 
gen in the inspired gas for breath n during the first twenty 
breaths the subjects are breathing air, (e.g., inspired 
F,O,=0.21), and are then changed to a gas mixture enriched 
in oxygen (e.g., F,O,=1.0). 

At the end of the first breath following the switch in F,O,, 
the concentration in the unit is: 


Vo-Co + AV-C! (5) 


insp 


i Vo +AV 


The same reasoning can be used to establish a recursive 
formula for concentration after n breaths: 


Vo Cyt +AV-C? (6) 


insp 


a Vo + AV 


which in turn can be re-written as a function of the unit’s 
specific ventilation, with SV=AV/V,, as: 


(7) 


= Gi 


Eq. (7) can be used to model oxygen concentration. A 
simulated response for an individual voxel can be obtained 
for a given set of inspired O, fraction breaths (e.g., inspired 
oxygen fraction (650), as seen in the dotted line in FIG. 6). 
For example, FIG. 6 also presents the time series of two 
units: unit (630) with low specific ventilation (e.g., SV=0.2) 
and unit (640) with high specific ventilation (e.g., SV=0.8). 

A linear relationship can exist between R=1/T, and F,O,. 
The slope and zero crossing values reported for the R1-F,O, 
function can result in a almost linear relationship between T, 
and F,O, in the range used in the exemplary model and 
simulation (e.g., F}O,=0.21 or 1.0). 

Based on numerical simulations, for the inversion time 
described (e.g., 1000 msec) and the variation of T1 with 
inspired oxygen concentration, the error involved in assum- 
ing that the MR signal varies linearly with oxygen concen- 
tration is <3% for these exemplary implementations. The 
calibration curve can be based on this assumption. 

Time-Shifted Cross Correlation Between the Driving 
Function and Simulated Units 

To determine the correlation delay time, the cross corre- 
lation between each simulated ventilatory unit response and 
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the inspired oxygen fraction can be computed and is repre- 
sented in FIG. 6 (e.g., simulated ventilatory unit response for 
unit (630) and unit (640) and inspired oxygen fraction 
(650)). The inspired F,O, curve can be shifted breath-by- 


breath, and the time shift that maximizes the time-shifted 5 


crossed correlation can be a measure of how fast the unit 
equilibrates, e.g., correlation delay time. 

FIG. 7A shows exemplary plot (701) featuring correlation 
delay time. Correlation delay time can be referred to as the 
time delay that maximizes the cross-correlation of the mea- 
sured signal time course from an image voxel with a square 
wave representing the onset of increased inspired O,. Plot 
(701) shows two exemplary correlation delay times (corre- 
lation delay time (702) for 2 breaths and correlation delay 
time (703) for 6 breaths) determined based on measured 
signal time course (704) and rise time (705). Correlation 
delay time can be used as an empirical measure of the rise 
time (e.g., rise time (705) shown as dashed lines in FIG. 7A). 
For example, mathematically, the cross correlation is maxi- 
mized when the response is shifted by about half the rise 
time, so that different rise times effectively translate to 
different delays in the cross-correlation analysis. This allows 
for the use of correlation delay time as an empirical index of 
the rise time. This calculation can be done for the entire time 
series corresponding to each image voxel. It is noted that in 
some exemplary systems, a dead space, for example, in the 
plumbing leading from the gas containers to the subject, can 
introduce a true global delay, but this delay can be calculated 
from the geometry and flow rate of the delivery system, and 
eliminated before the cross-correlation analysis. 


Age Ht. Wt 
(yrs) 


Subj Gndr 
S1 M 33 
$2 M 26 
$3 F 26 
$4 F 39 
$5 F 24 
$6 M 52 
$7 M 34 
$8 M 28 1.70 60 


Mean+SD 33+9 


1.76+0.08 82217 3.85 + 0.61 
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and the simulated signal was analyzed using the time-shifted 
cross correlation. The outcome of these exemplary simula- 
tions for units with specific ventilation in the range specified 
is shown in plot (751) of FIG. 7B. Plot (751) can include a 
conversion tool that allows the translation of correlation 
delay time (e.g., measured using the proton-MRI acquired 
time series) into quantitative values of ventilation, on a 
voxel-per-voxel basis. 

Plot (751) shows that for a given measured correlation 
delay time (shown on the y-axis in FIG. 7B), the corre- 
sponding specific ventilation can depend not only on the 
specific ventilation, but also on an extrinsic delay (for 
example, which can result from plumbing volume in the 
inspiratory line, discussed below). In essence, any volume 
that exists within the experimental configuration between a 
subject’s mouth and the valve used to change between gases 
of different F,O, can introduce a delay that can appear to 
artificially reduce SV unless, properly accounted for. Thus, 
it is beneficial to determine this delay by measuring respi- 
ratory flow and inspiratory plumbing volume. 

Exemplary MRI Data Collection Subjects 

In exemplary experiments implementing the exemplary 
method to quantify specific ventilation (SV), eight healthy 
subjects (e.g., 3 female, 5 male subjects) were studied in the 
supine posture. Table 1 presents subject characteristics (e.g., 
gender, age, height and weight) and pulmonary function data 
(e.g., FEV,, FVC, and FEV,/FVC). Table 1 shows exem- 
plary subject characteristics, pulmonary function data and 
slope of the lung height—Specific ventilation relationship 
for the eight exemplary subjects. 


TABLE 1 

FEV1 FVC 
(liters) (liters) FEV1/FVC Spatial 

(% (% (% -1/slope fractal 

predicted) predicted) predicted) (em!) dimension 

3.99 5.17 0.77 0.020 1.15 
88% 93% 95% 

4.57 5.28 0.87 0.015 1.15 
92% 87% 05% 

4.13 5.04 0.82 0.031 15 
115% 118% 96% 

3.03 3.65 0.83 0.038 09 
101% 99% 01% 

3.08 3.65 0.84 0.043 stl 
102% 99% 02% 

4.50 5.36 0.84 0.046 10 
105% 96% 09% 

3.84 4.98 0.77 0.021 16 
89% 93% 95% 

4.43 5.16 0.86 0.019 3 
107% 103% 04% 

4.79 £0.71 0.82 +£0.04 0.029* + 0.012 1.13 + 0.03 
100 + 10% 99 + 9% 101 + 5% 


*Group average significantly different from zero, p = 0.0002. 


To relate the correlation delay time to the specific venti- 55 


lation, an exemplary model of a lung unit can be imple- 
mented that simulates the series of consecutive breaths 
following a sudden change in inspired F,O,. The simulated 
data for units with different specific ventilations can be used 


to calibrate an exemplary signal-processing algorithm, in 60 


which the rapidity of the increase/decrease in the MRI signal 
is measured as a functional delay between the change in 
F,O, and the resulting time course of the MRI signal 
intensity. The measured functional delay can be translated 


into quantitative specific ventilation. 65 


Exemplary simulations were performed in units with 
specific ventilation ranging from 0.05 to 1 in steps of 0.05, 


MRI Materials and Equipments 

In exemplary experiments implementing the exemplary 
method to quantify specific ventilation (SV), MRI data was 
collected from the exemplary subjects using a 1.5 Tesla 
Signa HDx TwinSpeed MRI system (e.g., General Electric 
Medical Systems, Milwaukee, Wis., USA). A single sagittal 
slice was selected in the right lung, e.g., in order to avoid 
physiological noise arising from cardiac movements. Slice 
selection was aimed at selecting the slice within the lung 
presenting the largest anterior-posterior dimension, while 
avoiding major hilar vessels. 

Exemplary Specific Ventilation Imaging (SVI) Proce- 
dures 
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In exemplary experiments implementing the exemplary 
method to quantify specific ventilation (SV), two-dimen- 
sional T, weighted images were acquired using an inversion 
recovery (T,=1000 msec) single shot fast spin echo (SSFSE) 
sequence, with images being acquired with a half Fourier 
acquisition (HASTE), and with a 40x40 cm field of view, 
echo time of approximately 30 ms, and a 15 mm image slice 
thickness. A MRI homodyne reconstruction algorithm was 
used to rescale the data to a 256x256 matrix, with each voxel 
thus corresponding to ~1.6x1.6x15 mm (~40 mm*). An 
inversion time (T,) of 1000 msec, approximating T, of lung 
tissue, was implemented to ensure maximal sensitivity to 
changes in O, concentration. The long repetition time used 
in this single shot fast spin echo sequence (e.g., 5 s, 
compared to the T, of blood, 1.4 s) can render the O, 
induced contrast independent of lung density. However, lung 
density does alter the local signal-to-noise ratio. In the 
exemplary implementations, HASTE acquisition was uti- 
lized to keep the echo time short to minimize signal loss 
observed in the lung. Images were voluntarily respiratory 
gated. For example, subjects were instructed to take a 
normal breath in following the noise made by the MR image 
acquisition, and relax back to Functional Residual Capacity 
(FRC) at a comfortable expiratory flow rate. Subjects were 
comfortable with the default 5 s inter-breath intervals (12 
breaths per minute). MR images were acquired during a 
short (e.g., ~hundreds of milliseconds), post-expiratory 
breath hold at FRC resulting in a relatively natural respira- 
tory maneuver, which was only constrained by a constant 
breathing rate. 

In exemplary experiments implementing the exemplary 
method to quantify specific ventilation (SV), the addition of 
O, as a contrast agent followed a block design functional 
Magnetic Resonance Imaging (f{MRI) approach, which can 
be utilized for neuroimaging studies. An MR image ofa lung 
was acquired every 5 s, e.g., with 20 images acquired with 
the subject inspiring air (21% oxygen) and 20 images while 
inspiring 100% O2. This cycle was repeated 5 times, and an 
additional 20 breaths of 100% oxygen were added at the end 
of the last cycle making a total of 220 images (total imaging 
time was 18 min 20 sec) for an exemplary experiment 
involving each subject. Blocks of 20 breaths of air and 100% 
oxygen can ensure an approach to full equilibration for the 
ranges of specific ventilation and V_,/Q ratios observed in 
normal healthy subjects. Five cycles of air-oxygen can 
provide an acceptable signal to noise ratio. While more or 
less cycles of air-oxygen can be implemented, the exemplary 
five cycles of air-oxygen were implemented as a compro- 
mise between keeping the total acquisition time below 20 
min and an improved signal to noise ratio that would result 
from a longer sequence. 

In exemplary experiments implementing the exemplary 
method to quantify specific ventilation (SV), a face mask 
(e.g., Hans Rudolph, Kans., USA, with deadspace of 73 to 
113 mL, depending on mask size) equipped with a non- 
rebreathing T-valve (e.g., with deadspace of 27.9 mL) was 
fitted to the subject. One end of the T-valve was connected 
to the inlet, where a remote controlled 3-way pneumatic 
sliding valve (e.g., Hans Rudolph, model 8500) allowed 
rapid switching between room air and O, contained in a 170 
liter Douglas bag (e.g., Hans Rudolph type 6170). The 
inspiratory path resistance on the room air and O, circuits 
were matched to eliminate changes in FRC following 
changes in inspiratory path. The outlet of the T-valve was 
connected to a ~6 m long large bore low resistance expira- 
tory line, leading out of the scanner room, where expired 
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tidal volume was simultaneously measured using a Parvo- 
Medics Metabolic Measurement System (ParvoMedics, 
Sandy, Utah, USA). 

The exemplary experimental configuration for controlling 
the inspired gas (room air, 100% oxygen) introduced an 
extrinsic plumbing delay in the inspired signal that can be 
determined by the volume of tubing that connects the remote 
control valve to the subject. This volume was made as small 
as possible (e.g., 0.6 L) under constraints imposed by 
working in an MRI environment. Neglecting diffusion, this 
delay can be computed from each subject’s tidal volume and 
the tubing volume, by dividing the tubing volume by the 
tidal volume. For example, in practice, this extrinsic delay 
can be taken into account by computing individual inspired 
fractional oxygen concentration (F,O,) time series on a 
breath-by-breath basis (e.g., Cries in Eq. (7)), taking the 
tubing volume as a delay chamber of fixed volume, through 
which the inspired gas must pass. Once this extrinsic plumb- 
ing delay is corrected for, what is left corresponds to the 
signal intensity change over time, allowing the computation 
of specific ventilation as described in the next section. 

Exemplary Specific Ventilation Imaging (SVI) Data 
Analysis 

In exemplary experiments implementing the exemplary 
method to quantify specific ventilation (SV), a time series of 
the exemplary 220 images corresponding to the SVI 
sequence was constructed. Quality control of the acquired 
image can be performed at this stage. For example, images 
in which a subject was not at FRC based on diaphragm 
position compared to adjacent images can be removed from 
the series and replaced by an interpolated image, constructed 
from the preceding and following images. In exemplary 
experiments, a region of interest was manually drawn, and 
the subsequent analyses were restricted to the voxels inside 
the lung. For example, the region of interest encompassed 
the entire lung, yet avoided partial volume effects from 
regions close to the chest wall and the diaphragm. Data 
analysis was performed on a voxel-by-voxel basis, the time 
course of each voxel (e.g., MRI signal intensity versus time) 
was used to compute the regional correlation delay time. For 
example, data analysis can be performed using Matlab 
software (Mathworks, Natick, Mass.). 

In exemplary experiments implementing the exemplary 
method to quantify specific ventilation (SV), the correlation 
delay time was computed using a shifted cross correlation 
between the inspired F,O, and the acquired signal intensity 
time series. The time shifted cross correlation can be con- 
sidered as a fast and simple approach to implement, fitting 
only one parameter (e.g., the correlation delay time), and is 
independent of the asymptotic signal intensity. For example, 
by removing this extra degree of freedom, the model can 
remain simple while capturing only the time course of the 
transition. This can leave out, for example, the physiological 
dead space, water vapor, ratios; factors that alter the steady 
state CO, concentration, and any effect of varying V_,/Q 
concentration, but not the time course of the transition. 
Moreover, the shifted correlation delay approach can act as 
a smoothing low pass filter, eliminating some of the high 
frequency noise present in the exemplary data that did not 
allow a direct fit of the model (Eq. (7)). In practice, for 
example, the inspired F,O, (driving function) can be corre- 
lated with the time course of each voxel, and this process can 
be repeated for delayed versions of the driving function 
(e.g., delayed by an integer or fractional number of breaths). 
The delay that maximizes the cross correlation between the 
time shifted driving function and the actual voxel response 
is the correlation delay time. The correlation delay time 
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computed in this exemplary way can retain voxels whose 
correlation with the optimal shifted inspired fractional oxy- 
gen concentration is significant (P<0.05), and the null 
hypothesis of no correlation can be rejected in these cases. 
In voxels in which the null hypothesis is accepted, the 
corresponding lung voxel can be attributed no specific 
ventilation value and treated as missing data. 

FIG. 8 shows an exemplary time course plot (800) of 
signal intensity (801), represented by filled circles, continu- 
ous lines) for one such exemplary voxel. The arbitrarily 
chosen exemplary voxel plotted was located in the depen- 
dent portion of the lung. The corresponding driving function 
(802) representing F,O, is also shown, represented by the 
dotted line. The delay between the driving function and the 
measured signal intensity is a measure of the correlation 
delay time for that voxel. 

Exemplary Statistical Analyses 

In exemplary experiments implementing the exemplary 
method to quantify specific ventilation (SV), each series of 
220 breaths was considered as a single measure of specific 
ventilation, computed as described above, creating for each 
subject one map of specific ventilation. All the voxels 
presenting a statistically significant (p<0.05) correlation 
with the optimally shifted drive function were included in 
the exemplary analysis. The remaining voxels were treated 
as missing data. Two different exemplary analyses were 
implemented, e.g., a comparison of specific ventilation by 
lung thirds and a finer analysis on a cm-per-cm basis moving 
up the supine lung. 

For the exemplary lung thirds analysis, specific ventila- 
tion was partitioned into three gravitational regions, corre- 
sponding to thirds of the lung based on equal vertical extent, 
e.g., the dependent portion, the intermediate region, and the 
non-dependent region. The data was reduced to a subject- 
by-subject average specific ventilation per lung gravitational 
region. One-way ANOVA was used to compare the gravi- 
tational gradient in specific ventilation across the lung 
regions (e.g., the 3 levels: dependent, intermediate, and 
non-dependent). Where overall significance was present, 
post-hoc testing was conducted using Student’s t-test. 

For the exemplary finer analysis (cm-per-cm basis mov- 
ing up the supine lung), linear regression was used to 
evaluate the linear relationship between the vertical height 
of the lung as a function of specific ventilation, e.g., by 
dividing the data into isogravitational slices of 1 cm thick- 
ness (~6 vertical voxels). As the exemplary subjects were 
studied in the supine position, the vertical height of the lung 
(isogravitational level) was measured along the anterior- 
posterior axis, with zero height corresponding to the most 
dependent isogravitational voxels. The linear relationships 
were evaluated individually for each subject. The inter- 
subject averaged specific ventilation versus lung height 
values correspond to averages over vertical 1 cm regions for 
the 8 subjects. 

It is noted that different subjects have different anterio- 
posterior lung dimensions; therefore, results were reported 
only for lung heights that included data from the eight 
subjects. The slope of the individual relationships between 
height and specific ventilation were compared to a zero 
slope using a one-group t-test. 

It is noted that exemplary data are presented as meanzSD. 
In the case in which data for the eight subjects were 
averaged, SD corresponds to the inter-subject variability 
(SD). Regarding the calculating relative dispersion, spatial 
SD refers to the inter-voxel variability within a specific 
ventilation map. The null hypothesis (no effect) was rejected 
when P<0.05, two tailed, except where otherwise indicated. 
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Exemplary statistical analyses were performed using Prism 
(GraphPad, San Diego, Calif.). 

Exemplary Results 

General results of the exemplary experiments implement- 
ing the exemplary method to quantify specific ventilation 
(SV) included the following. Subject descriptive data and 
pulmonary function measurements are presented in Table 1. 
The subjects studied in the exemplary experiments had 
normal spirometry, as indicated by an average 
FEV ,=100+10%-predicted, FVWC=99+9%-predicted, and 
FEV ,/FVC=101+5%-predicted (as shown in Table 1). The 
exemplary subjects maintained a relatively constant expired 
tidal volume throughout the experiment, e.g., averaging 
0.87+0.18 L. Heart rate during SVI data acquisition aver- 
aged 65+9 beats/min (e.g., no difference between inspiring 
air and 100% oxygen), and arterial oxygen saturation mea- 
sured by pulse oximetry was 97.4+1.7% while breathing 
room air and increased to 97.8+1.5% during inspiration of 
100% oxygen. The selected breathing rate was followed by 
all subjects in the exemplary experiments, and tidal volume 
was relatively stable throughout the 220 breaths (e.g., the 
average CV for the eight subjects was 16%). Tidal volume 
did not change when comparing breaths taken while inspir- 
ing air and oxygen (e.g., average over the eight exemplary 
subjects, 0.86+0.15 Liters on air, 0.87+0.17 Liters, 2-tailed 
t-test paired comparison, p=0.82). 

SVI results of the exemplary experiments implementing 
the exemplary method to quantify specific ventilation (SV) 
included the following. FIG. 8 presents the time series (e.g., 
time course plot (800)) of the signal intensity of a voxel 
(e.g., signal intensity (801)) together with the inspired F,O, 
(function (802)) for a voxel in the dependent lung region. 
The exemplary experiment included five repetitions of a 
block of 20 air breaths and 20 100%-oxygen breaths per- 
formed by a subject. As shown in the figure, when the 
inspired line was changed from air to 100% oxygen, signal 
intensity increased; signal intensity decreased for the oppo- 
site change, e.g., from 100% oxygen back to room air. After 
the 5 cycles of air and 100% oxygen, twenty additional 
breaths on 100% oxygen were added in the exemplary 
experiment. Noise can be expressed as the coefficient of 
variation of the steady state signal (e.g., first 20 breaths on 
air, and last 20 breaths inspiring O.) averaged over the entire 
lung of the 8 subjects studies. The exemplary noise was 
determined as 22.841.5% for air, and 18.4+1.8% for oxygen. 
In order to test the robustness of the exemplary algorithm, 
the disclosed model was run with similar levels of random 
noise. The average specific ventilation estimation error 
introduced by 20% noise levels can be considered as 
extremely small (e.g., average error in estimating SV 
~0.003). The maximal SV estimation error was determined 
to be ~0.013 and was observed for very low specific 
ventilation units (SV<0.02), e.g., which is below the normal 
healthy range. 

From each individual voxel response, a map of the 
correlation delay time for all voxels within the lung can be 
computed. An exemplary individual voxel response maps is 
presented for a representative subject, as shown by map 
(901) in FIG. 9A. The arrow shown indicates the direction 
of gravity; the head is located to the right of the image, and 
the diaphragm to the left. In FIG. 9A, warmer colors (e.g., 
colors trending toward red) represent portions of the lung 
with a longer correlation delay time. As previously dis- 
cussed, the extrinsic delay resulting from the inspiratory 
plumbing was accounted for prior to the analysis. 

Using the results presented in FIG. 7B, this correlation 
delay time was then translated into a quantitative measure of 
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specific ventilation, as shown by map (902) in FIG. 9B. In 
FIG. 9B, warmer colors (e.g., colors trending toward red) 
represent regions of the lung with higher specific ventilation. 
The quantitative measure of specific ventilation results of 
the exemplary experiments implementing the disclosed 
technology include average specific ventilation in a slice of 
the right lung (averaged over all subjects) determined to be 
0.33+0.11. Overall specific ventilation heterogeneity, as 
measured by the relative dispersion, averaged 0.63+0.11; 
with an individual range of 0.50 to 0.79. 

The exemplary results of quantitative measure of specific 
ventilation as shown in FIG. 9B demonstrates a clear vertical 
(gravitational) gradient of specific ventilation is present. 
Units in the dependent portions of the lung have higher 
specific ventilation than those in the non-dependent portions 
of the lung. This relationship is quantified in graph (1000) of 
FIG. 10, which is shown in the results for the regional 
dependence of specific ventilation over the eight exemplary 
subjects in the three regions, e.g., nondependent region 
(1001), intermediate region (1002) and dependent region 
(1003). For example, a statistically significant difference in 
specific ventilation was observed, with the most dependent 
third of the lung having a specific ventilation of 0.42+0.14, 
the intermediate third 0.29+0.10 and the non-dependent 
third 0.24+0.08 (all differences, P<0.05). 

An approximately linear decrease in specific ventilation 
with height was observed, as evidenced in FIGS. 11A-11C. 
FIG. 11A shows specific ventilation plotted against height 
for one exemplary subject in plot (1110). FIG. 11B shows 
specific ventilation plotted against height for each of the 
eight exemplary individual subjects in plot (1120). FIG. 11C 
shows the average specific ventilation plotted against height 
of all eight exemplary subjects in plot (1130). The inverse of 
the slope of the plotted line shown in FIGS. 11A-11C is a 
measure of the decrease in specific ventilation per cm. 
Specific ventilation was found to decrease 0.029+0.012/cm 
on average for the eight exemplary individual subjects, 
determined from the inverse of the (slope 1131) shown in 
FIG. 11C. 

The exemplary implementation of the disclosed O,-en- 
hanced proton MR imaging technology showed regional 
quantification of specific ventilation, which can, for 
example, extract more information than just absence or 
presence of ventilation. For example, from a physiological 
standpoint, specific ventilation imaging (SVI) can show a 
clear vertical (gravitational dependent) gradient with higher 
values, which was demonstrated in the exemplary experi- 
ments as shown in dependent lung in FIG. 10. 

The spring like self-deformation of the lung, e.g., the 
Slinky effect, is likely the main determinant of the gravita- 
tional gradient observed in specific ventilation. However, in 
the most dependent third of the supine lung there are 
deviations from this overall behavior suggesting important 
influences of other factors such as the non-linearity of the 
pressure-volume curve of the lung, the large vessels in the 
lung, and dependent airways closure. 

The exemplary measurements show a linear vertical gra- 
dient in specific ventilation (FIGS. 9A-4C and 11A and 
11B). Dependent portions of the lung have a higher specific 
ventilation that the non-dependent portions (FIG. 10). The 
dependent portion of the lung partially supports the weight 
of the upper portions, and thus is more “compressed.” 
Therefore its end-expiratory volume is smaller. Further, the 
dependent portion is subject to a higher (less negative) 
pleural pressure than an equivalent non-dependent region, 
therefore placing it on a steeper portion of the pressure- 
volume curve. For example, a change in pleural pressure 
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resulting from an inspiratory effort may result in a greater 
increase in volume. Both the lower initial (end-expiratory) 
volume and the increased change in volume can be expected 
to contribute to the higher specific ventilation observed in 
the dependent regions of the lung, when compared to the 
non-dependent. 

Quantification of specific ventilation in the human lung 
can be implemented by the described method of proton-MRI 
SVI, which uses O, as an inspired contrast agent. 

Implementations of the subject matter and the functional 
operations described in this specification can be imple- 
mented in various systems, digital electronic circuitry, or in 
computer software, firmware, or hardware, including the 
structures disclosed in this specification and their structural 
equivalents, or in combinations of one or more of them. 
Implementations of the subject matter described in this 
specification can be implemented as one or more computer 
program products, i.e., one or more modules of computer 
program instructions encoded on a tangible and non-transi- 
tory computer readable medium for execution by, or to 
control the operation of, data processing apparatus. The 
computer readable medium can be a machine-readable stor- 
age device, a machine-readable storage substrate, a memory 
device, a composition of matter effecting a machine-read- 
able propagated signal, or a combination of one or more of 
them. The term “data processing apparatus” encompasses all 
apparatus, devices, and machines for processing data, 
including by way of example a programmable processor, a 
computer, or multiple processors or computers. The appa- 
ratus can include, in addition to hardware, code that creates 
an execution environment for the computer program in 
question, e.g., code that constitutes processor firmware, a 
protocol stack, a database management system, an operating 
system, or a combination of one or more of them. 

A computer program (also known as a program, software, 
software application, script, or code) can be written in any 
form of programming language, including compiled or 
interpreted languages, and it can be deployed in any form, 
including as a stand alone program or as a module, compo- 
nent, subroutine, or other unit suitable for use in a computing 
environment. A computer program does not necessarily 
correspond to a file in a file system. A program can be stored 
in a portion of a file that holds other programs or data (e.g., 
one or more scripts stored in a markup language document), 
in a single file dedicated to the program in question, or in 
multiple coordinated files (e.g., files that store one or more 
modules, sub programs, or portions of code). A computer 
program can be deployed to be executed on one computer or 
on multiple computers that are located at one site or dis- 
tributed across multiple sites and interconnected by a com- 
munication network. 

The processes and logic flows described in this specifi- 
cation can be performed by one or more programmable 
processors executing one or more computer programs to 
perform functions by operating on input data and generating 
output. The processes and logic flows can also be performed 
by, and apparatus can also be implemented as, special 
purpose logic circuitry, e.g., an FPGA (field programmable 
gate array) or an ASIC (application specific integrated 
circuit). 

Processors suitable for the execution of a computer pro- 
gram include, by way of example, both general and special 
purpose microprocessors, and any one or more processors of 
any kind of digital computer. Generally, a processor will 
receive instructions and data from a read only memory or a 
random access memory or both. The essential elements of a 
computer are a processor for performing instructions and 
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one or more memory devices for storing instructions and 
data. Generally, a computer will also include, or be opera- 
tively coupled to receive data from or transfer data to, or 
both, one or more mass storage devices for storing data, e.g., 
magnetic, magneto optical disks, or optical disks. However, 
a computer need not have such devices. Computer readable 
media suitable for storing computer program instructions 
and data include all forms of non volatile memory, media 
and memory devices, including by way of example semi- 
conductor memory devices, e.g., EPROM, EEPROM, and 
flash memory devices. The processor and the memory can be 
supplemented by, or incorporated in, special purpose logic 
circuitry. 

While this patent document contains many specifics, these 
should not be construed as limitations on the scope of any 
invention or of what may be claimed, but rather as descrip- 
tions of features that may be specific to particular embodi- 
ments of particular inventions. Certain features that are 
described in this patent document in the context of separate 
embodiments can also be implemented in combination in a 
single embodiment. Conversely, various features that are 
described in the context of a single embodiment can also be 
implemented in multiple embodiments separately or in any 
suitable subcombination. Moreover, although features may 
be described above as acting in certain combinations and 
even initially claimed as such, one or more features from a 
claimed combination can in some cases be excised from the 
combination, and the claimed combination may be directed 
to a subcombination or variation of a subcombination. 

Similarly, while operations are depicted in the drawings in 
a particular order, this should not be understood as requiring 
that such operations be performed in the particular order 
shown or in sequential order, or that all illustrated operations 
be performed, to achieve desirable results. Moreover, the 
separation of various system components in the embodi- 
ments described above should not be understood as requir- 
ing such separation in all embodiments. 

Only a few implementations and examples are described 
and other implementations, enhancements and variations 
can be made based on what is described and illustrated in 
this patent document. 

What is claimed is: 

1. Acomputer-implemented method using magnetic reso- 
nance imaging to characterize ventilation and perfusion in a 
lung, comprising: 

acquiring, using a magnetic resonance imaging (MM) 

machine under the control of a processing unit, a 
magnetic resonance (MR) image of the lung that 
includes MR data in a voxel of the MR image, wherein 
the acquiring causes the MM machine to measure a 
breathing frequency value; 

determining, by the processing unit a water density value 

and a perfusion value in at least one voxel of the MR 
image based on the MR data; 
determining, by the processing unit from the water den- 
sity value, an air content value in the at least one voxel; 

determining, by the processing unit, a specific ventilation 
value in the at least one voxel based on a time delay 
between an onset of a stimulus to the lung and a 
response detected in the at least one voxel; 

determining, by the processing unit a ventilation-perfu- 
sion ratio value that is the product of the specific 
ventilation value, the air content value, the inverse of 
the perfusion value, and the breathing frequency value; 
and 

evaluating, by the processing unit, a spatial distribution of 

the specific ventilation value, perfusion value, and the 
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ventilation-perfusion ratio by computing a standard 

deviation of the specific ventilation value and perfusion 

value; 

wherein the specific ventilation value is computed by 
SVI*(1-density)*f, wherein SVI represents quanti- 
tative data from the MR image, (1-density) repre- 
sents the air content value, and f, represents the 
breathing frequency value, thereby determining 
when the lung is diseased. 

2. The method of claim 1, wherein the MR image is a 2D 
image or a 3D image. 

3. The method of claim 1, wherein the determining the 
ventilation-perfusion ratio value is determined voxel-by- 
voxel. 

4. The method of claim 1, wherein the acquiring the MR 
image uses oxygen as an inhaled contrast agent. 

5. The method of claim 1, further comprising determining 
an alveolar ventilation value that is the product of the 
specific ventilation value, the air content value, and the 
breathing frequency value. 

6. The method of claim 5, further comprising: 

computing a first variance in the specific ventilation 

value, a second variance in the perfusion value, and a 
third variance in the ventilation-perfusion ratio value; 
and 

generating a correlation coefficient from one or more of 

the first variance, the second variance, and the third 
variance. 

7. The method of claim 1, further comprising producing 
an image that shows the spatial distribution of ventilation, 
perfusion, and ventilation-perfusion for each voxel associ- 
ated with the acquired MR image. 

8. The method of claim 1, wherein the spatial distribution 
of ventilation, perfusion, and ventilation-perfusion is used as 
an indicator of a diseased lung. 

9. The method of claim 1, further comprising: 

determining ventilation-perfusion matching in the lung 

from at least the ventilation-perfusion ratio value. 

10. The method of claim 9, wherein ventilation-perfusion 
matching is used as an indicator of a diseased lung. 

11. A magnetic resonance imaging system to characterize 
ventilation and perfusion in a lung, comprising: 

a magnetic resonance imaging machine that acquires an 

MR image of the lung; and 

a processing unit having a processor and a non-transitory 

computer-readable storage medium having instructions 

stored thereon that when executed by the processor 

cause the processing unit to at least: 

command the magnetic resonance imaging machine to 
acquire the MR image that includes MR data in each 
voxel of the MR image, 

process the MR data to determine a water density value 
and a perfusion value in at least one voxel of the MR 
image, wherein the water density value is used to 
determine an air content value, 

measure a breathing frequency of the lung; 

determine a specific ventilation value in the at least one 
voxel based on a time delay between an onset of a 
stimulus to the lung and a response detected in the at 
least one voxel; 

determine a ventilation-perfusion ratio value that is the 
product of the specific ventilation value, the air 
content value, the inverse of the perfusion value, and 
the breathing frequency; and 

evaluate a spatial distribution of the specific ventilation 
value, perfusion value, and the ventilation-perfusion 
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ratio by computing a standard deviation of the spe- 

cific ventilation value and perfusion value; 

wherein the specific ventilation value is computed by 
SVI*(1-density)*f, wherein SVI represents quan- 
titative data from the MR image, (1-density) 
represents the air content value, and f, represents 
the breathing frequency value, thereby determin- 
ing when the lung is diseased. 

12. The magnetic resonance imaging system of claim 11, 
wherein the MRI system uses oxygen as an inhaled contrast 
agent to acquire the MR image of the lung. 

13. The magnetic resonance imaging system of claim 11, 
wherein the instructions, when executed by the processor, 
further cause the processing unit to at least: 

determine an alveolar ventilation value that is the product 

of the specific ventilation value, the air content value, 
and the breathing frequency. 
14. The magnetic resonance imaging system of claim 13, 
wherein the instructions, when executed by the processor, 
further cause the processing unit to at least: 
compute a first variance in the specific ventilation value, 
a second variance in the perfusion value, and a third 
variance in the ventilation-perfusion ratio; and 

generating a correlation coefficient from one or more of 
the first variance, the second variance, and the third 
variance. 

15. The magnetic resonance imaging system of claim 11, 
wherein the instructions, when executed by the processor, 
further cause the processing unit to at least: 

produce an image that shows the spatial distribution of 

ventilation, perfusion, and ventilation-perfusion for 
each voxel associated with the acquired MR image. 

16. The magnetic resonance imaging system of claim 11, 
wherein the processing unit uses the spatial distribution of 
ventilation, perfusion, and ventilation-perfusion as an indi- 
cator of a diseased lung. 

17. The magnetic resonance imaging system of claim 11, 
wherein the instructions, when executed by the processor, 
further cause the processing unit to at least: 

determine ventilation-perfusion matching in the lung 

from at least the ventilation-perfusion ratio value. 
18. The magnetic resonance imaging system of claim 17, 
wherein the processing unit uses ventilation-perfusion 
matching as an indicator of a diseased lung. 
19. A computer program product comprising a non- 
transitory computer-readable storage medium having 
instructions stored thereon, the instructions comprising: 
code for acquiring, using a magnetic resonance imaging 
(MRI) machine, an MR image of a lung that includes 
MR data in a voxel of the MR image; 

code for determining a water density value and a perfu- 
sion value in at least one voxel of the MR image based 
on the MR data, wherein the water density value is used 
to determine an air content value; 

code for measuring a breathing frequency value; 


15 


30 


35 


40 


45 


28 


code for determining a specific ventilation value in the at 
least one voxel based on a time delay between an onset 
of a stimulus to the lung and a response detected in the 
at least one voxel; 

code for determining a ventilation-perfusion ratio value 

that is the product of the specific ventilation value, the 
air content value, the inverse of the perfusion value, 
and the breathing frequency value; and 

code for evaluating a spatial distribution of the specific 

ventilation value, perfusion value, and the ventilation- 

perfusion ratio by computing a standard deviation of 

the specific ventilation value and perfusion value, 

wherein the specific ventilation value is computed by 
SVI*(1-density)*f, wherein SVI represents quanti- 
tative data from the MR image, (1-density) repre- 
sents the air content value, and f, represents the 
breathing frequency value, thereby determining 
when the lung is diseased. 

20. The non-transitory computer-readable storage 
medium of claim 19, wherein the code for acquiring the MRI 
image of the lung includes code that uses oxygen as an 
inhaled contrast agent. 

21. The non-transitory computer-readable — storage 
medium of claim 19, further comprising code for determin- 
ing an alveolar ventilation value that is the product of the 
specific ventilation value, the air content value, and the 
breathing frequency. 

22. The non-transitory computer-readable — storage 
medium of claim 21, wherein code for determining the 
alveolar ventilation value includes code for computing a first 
variance in the specific ventilation value, a second variance 
in the perfusion value, and a third variance in the ventilation- 
perfusion ratio; and 

generating a correlation coefficient from one or more of 

the first variance, the second variance, and the third 
variance. 

23. The non-transitory computer-readable — storage 
medium of claim 19, wherein code for determining the 
alveolar ventilation value further comprises code for pro- 
ducing an image that shows the spatial distribution of 
ventilation, perfusion, and ventilation-perfusion for each 
voxel associated with the acquired MR image. 

24. The non-transitory computer-readable storage 
medium of claim 19, wherein code for determining the 
alveolar ventilation value further comprises code for using 
the spatial distribution of ventilation, perfusion, and venti- 
lation-perfusion to indicate a diseased lung. 

25. The non-transitory computer-readable storage 
medium of claim 19, further comprising code for determin- 
ing ventilation-perfusion matching in the lung from at least 
the ventilation-perfusion ratio value. 

26. The non-transitory computer-readable storage 
medium of claim 25, wherein code for determining venti- 
lation-perfusion matching further comprises code for using 
ventilation-perfusion matching to indicate a diseased lung. 
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